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' A characteristic feature of the Portevin - Le ChateUer effect or the jerky flow is the stick-shp 

I nature of stress-strain curves which is believed to result from the negative strain rate dependence 

' of the flow stress. The latter is assumed to result from the competition of a few relevant time scales 

controlling the dynamics of jerky flow. We address the issue of time scales and its connection to 
the negative strain rate sensitivity of the flow stress within the framework of a model for the jerky 
flow which is known to reproduce several experimentally observed features including the negative 
strain rate sensitivity of the flow stress. We attempt to understand the above issues by analyzing 
the geometry of the slow manifold underlying the relaxational oscillations in the model. We show 
that the nature of the relaxational oscillations is a result of the atypical bent geometry of the slow 
I manifold. The analysis of the slow manifold structure helps us to understand the time scales op- 

erating in different regions of the slow manifold. Using this information we are able to establish 
connection with the strain rate sensitivity of the flow stress. The analysis also helps us to provide 
Ph I a proper dynamical interpretation for the negative branch of the strain rate sensitivity. 

PACS number(s): 83.50.-v, 83.50.By, 81, 05.45.-a, 62.20.-x 
> ■ I. INTRODUCTION 

oo ■ 

The Portevin - Le Chatelier effect Q] or the jerky flow has been an object of continued interest in materials science 
I for quite some time. The phenomenon refers to an instability seen in the form of repeated stress drops followed by 
periods of reloading observed when tensile specimens are deformed in a certain range of strain rates and temperatures 
1^. The effect is seen in maiiy interstitial and substitutional metallic alloys (commercial aluminium, brass, alloys of 
' aluminium and magnesium H ) . Each of the load drops is related to the formation and propagation of dislocation 
(— I : bands . The traditional picture of the instability is that it stems from dynamic interaction of mobile dislocations 
• '-H ' with solute atoms and is called as dynamic strain ageing . It is this that is expected to lead to negative strain rate 
^ ; sensitivity (SRS) of the flow stress 

>'...... 

. ^ , Plastic flow is intrinsically nonlinear and therefore methods of nonlinear dynamics have a natural role to play in 

■ understanding plastic instabilities • Use of these new techniques have led to insights which were hitherto not 

5h ] possible. The first attempt to look at the phenomenon from a nonlinear dynamical angle was taken by Ananthakrishna 
. . . 1 and coworkers |^ , which offers a natural basis for the description of the time dependent aspects of the PLC effect 
which were ignored in the earlier theories [^-0. Their model allows for explicit inclusion and interplay of different 
time scales inherent in the dynamics of dislocations. These authors show that the occurrence of the instability is a 
consequence of Hopf bifurcation as a function of the applied strain rate. Many known features of the PLC effect such 
as the existence of a window of strain rates and temperatures within which it occurs, etc., were correctly reproduced. 
Most importantly, and for the first time, the negative SRS was shown to emerge naturally in the model, as a result of 
nonlinear interaction of the participating defects. It also predicts the existence of chaotic stress drops in a range of 
strain rates, which has been recently verified p^-p^. Even the number of degrees of freedom estimated turn out to 
be the same as in the model offering justification for ignoring spatial degrees of freedom (See Ref. 16 also). Further 
dynamical analysis of the model for the creep case has shown that the temperature dependence of the strain bursts 
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is consistent with experimental findings |17 1 



The study of the PLC effect from a dynamical angle has been useful in elucidating several features; but it has also 
brought certain other issues into sharp focus which were hitherto not investigated in depth. This paper is intended 
to address one such issue related to the time scales relevant to the dynamics of the PLC effect within the framework 
of the above model. This is reflected in the two well known attributes of the PLC effect, namely, the negative strain 
rate behavior of the flow stress and the stick-slip or relaxational nature of the dynamics reflected in the stress time 
series. In order to motivate, we will present arguments showing that conflicting conclusions can be arrived at when 
one analyses this question starting from these two angles. 

We start with a discussion of the well accepted physical picture of the PLC effect namely the dynamic strain ageing. 
At a qualitative level, theories of strain ageing already have an implicit suggestion that the occurrence of the negative 
SRS is related to the competition of diffusive time scale and the waiting time of dislocations at obstacles |^,^, even 
though, there is no dynamics involved in these theories. The physical picture of strain ageing is as follows. At small 
velocities, solute atoms have enough time to diffuse to the temporarily arrested dislocations thus providing additional 
pinning thereby impeding their breakaway from localized obstacles. Due to the constant applied strain rate, the over- 
all stress to keep the dislocations moving increases bringing the stress to a threshold level beyond which dislocations 
break away [^U^. At high velocities, such additional pinning due to solute atoms is not possible since the waiting 
time of dislocations at obstacles is too short for the diffusion to occur. A schematic diagram of the SRS is shown in 
Fig. 1. In the language of the stick-slip dynamics, the branch A'B' corresponds to the stick-state and B'C to the 
slip-state (Fig. 1). The slope of stress verses velocity ('friction coefficient') at low velocities is much higher than that 
corresponding to high velocities since in the former case, solute atoms have to be dragged along with the dislocations, 
while in the latter case there is no solute atmosphere. Based on physical considerations, these two stable branches 
are assumed to be separated by an unstable branch with a negative slope to reflect the nonaccessible nature. 



The occurrence of the negative flow rate characteristic is not just limited to the PLC effect |119|,|20|]. With particular 
reference to the conceptual aspects of the negative branch, we cite two other mechanical systems, namely, the peeling 
of an adhesive tape and frictional sliding of a block of material over another |l^,|2^ which shows the inaccessible 
nature of the negative slope branch. However, in the case of the PLC effect, there has been attempts to obtain exper- 



imental points in this domain of strain rates 1 21 2^ which has lead to some confusion about the measurability of the 
negative slope branch of the SRS which will be discussed later (Section 5). Therefore, it is important to understand 
the meaning of the negative branch of SRS from a dynamical point of view with reference to the PLC effect which 
hopefully will lead to a better understanding of other stick-slip phenomenon. 

Returning to the PLC effect. Penning j23j was the first to recognize that the negative SRS could be used to explain 
the strain rate jumps observed in experiments. Subsequently, the negative SRS feature has been used as an input 
into several theories [^4|j25| ] . In particular, it has helped to successfully explain the nature of yield drops occurring in 
different regimes of strain rate and temperature Pertinent to the our discussion of time scales inherent to the 
PLC effect, we note that in such theories, two slow time scales corresponding to the two dissipative branches. A' B' 
and CD', show up along with two fast time scales corresponding to the jumps ( B'C and D' A') in the strain rates. A 
more direct reflection of the time scales inherent to the dynamics of the PLC effect can be deduced from stress-strain 
curve. 

To facilitate discussion of time scales involved in an experimental stress-time series, consider the so called machine 
equation written as 

O-a = K[ea - Cp] (1) 

where aa,ea, refer to the stress, applied strain rate and plastic strain rate respectively and k is the combined elastic 
constant of the machine and the sample. We note that only stress is monitored by the load sensing devise. Apart 
from this, it is possible to measure the plastic strain rate using strain gauges or using cinematographic techniques 
p6| . Using Eq.(l), we can now identify different time scales in an experimental curve. Consider a typical stress-strain 
curve for an applied strain rate of 8.3 x 10~^s~^ for the PLC effect in Cu — lOYoAl is shown in Fig. 2. (For our 
purpose, we will ignore the nonperiodic nature.) From the saw-tooth shape of the stress - strain series, two points 
emerge: (a) the positive slope of (Ja — ta curve is close to the elastic loading rate ( Kea), and (b) the duration of each 
stress drop is very short. From Eq.(l), we see that the stress drop duration is the time interval during which ep(t) 
larger than e^. We also note that the changes in slopes, when they occur, are abrupt (within the recording accuracy 
of 0.05s~^). Knowing that ep is proportional to the mobile dislocation density and using Eq. (1), we can see that 
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the mobile dislocation density should be nearly constant in the rising part oi a — Ea curve and therefore correspond 
to the stick-state. Further, the short duration of the stress drop should be a result of rapid multiplication of mobile 
dislocation and therefore corresponds to the slip-state. This must be followed by the process of immobilization of 
dislocations. However, the abrupt change in the slope ( from negative to positive) also implies that immobilization 
time scale is also fast. Indeed, as is clear from Fig. 2, it is not possible to separate out these two fast time scales. Thus, 
from the stress-strain curve, we see only one slow time scale and two fast time scales which is in apparent conflict 
with what was argued from the schematic diagram of the orbit A'B'C'D' in Fig. 1. This discussion raises several 
questions relating the origin of these time scales causing jumps in dislocations densities which needs to understood if 
the above inconsistency has to be resolved. Specifically: (i) What is dynamical mechanism which keeps the mobile 
dislocation density constant and in low levels for long intervals of time? (ii) What are the mechanisms for rapid 
multiplication and immobilization of mobile dislocations ? As we shall see, resolving these issues will also help us to 
interpret the negative SRS in an appropriate way. Further, associating various time scales with different branches of 
the SRS provides a better insight into the stick-slip dynamics of the PLC effect. 

Analysis of time scales can be best understood from a dynamical point of view. It is well known that relaxation oscil- 
lations are at the root of stick-slip behavior. One of the standard ways of understanding relaxation oscillations is by 
analyzing the slow manifold geometry p7|-^l|] of the underlying model. Following this, we shall attempt to understand 
the above issues from the point of view of relaxation oscillations. The paper is organized as follows. In section 2, we 
briefly introduce the model along with the known results. In section 3, we state some bifurcation features relevant for 
further discussion. In section 4, we show that the nature of relaxation oscillation in the model is atypical and is due 
to the atypical nature of bent nature of the slow manifold of the model. This analysis further helps us to understand 
the dynamical basis of different time scales relevant to the PLC effect. In section 5, we discuss the concept of negative 
SRS and its measurement in some detail to highlight the meaning of the negative branch. Using the geometry of the 
slow manifold, we calculate the dependence of stress on the plastic strain rate and show the connection between the 
various branches of the SRS and the time scales operating in different regions of slow manifold which in turn helps 
us to resolve the inconsistency of time scales. Section 6 is devoted to discussion and conclusions. 



II. A DYNAMICAL MODEL FOR JERKY FLOW 



The model consists of mobile dislocations and immobile dislocations and another type which mimics the Cottrell's 
type, which are dislocations with clouds of solute atoms Let the corresponding densities be Nm, Nim, and Ni, 
respectively. The rate equations for the densities of dislocations are: 

N„, = eVmN,n - PNl ~ l3N,rrN,r,^ + lN,„, - , (2) 

N^m = ml - m^mNm " 7^^m + a^N,, (3) 

N, = a„,N„, ~a,N,. (4) 

The overdot, here, refers to the time derivative. The first term in Eq. (2) is the rate of production of dislocations 
due to cross glide with a rate constant 9. Vm is the velocity of the mobile dislocations which in general depends 
on some power of the applied stress aa- The second term refers to the annihilation or immobilization of two mobile 
dislocations. The third term also represents the annihilation of a mobile dislocation with an immobile one. The fourth 
term represents the remobilization of the immobile dislocations due to stress or thermal activation ( see jNim in Eq. 
3). The last term represents the immobilization of mobile dislocations either due to solute atoms or due to other 
pinning centers. a„i refers to the concentration of the solute atoms which participate in slowing down the mobile 
dislocations. Once a mobile dislocation starts acquiring solute atoms we regard it as a new type of dislocation, namely 
the Cottrell's type Ni, i.e, the incoming term in Eq. (4). As they acquire more and more solute atoms they slow down 
and eventually stop the dislocation entirely. At this point, they are considered to have transformed to Nim (loss term 
in Eq. (4) and a gain term in Eq. (3)). Indeed, the whole process can be mathematically represented by defining 
Ni ~ K{t — t')N,n{t')dt' = am exp{—ai{t — t'))N„i{t')dt' which represents the entire process of slowing down 
of Nm in an exponential fashion with a time constant a^. (The choice of K as having a exponential form is obviously 
a simplification of the actual process.) 

These equations should be dynamically coupled to the machine equation which now takes the form 

da = nie'a - BoNmVm), (5) 
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where Vm is the velocity of mobile dislocations and Bq is the Burgers vector. A power law dependence of Vm = 
Vq (cTa/cTo)'" is used. These equations can be cast into a dimensionless form by using scaled variables x = Nm {P/l), 
V = Niyn [P/OVo), z = Ni {l3ai/jarn), t = 9Vot and = aa/a-Q. 

X — <j)^x — ax — b^x^ — xy + y, (6) 

y = bo (box'^ ~ xy - y + az) , (7) 

i = c{x — z), (8) 

^ = d (e - , (9) 

Here a = am/SVo,bo — j/6Vq,c — ai/9Vo, k — {9f3aod/jBo) and e = {e'aP/BoVoj). For these set of equations there 
is only one steady state which is stable. There is a range of the parameters a, 6, c, d, m and e for which the linearized 
equations are unstable. In this range x, y, z and are oscillatory. 



Among these physically relevant parameters, we study the behavior of the model as a function of most important 
parameters namely the applied strain rate e and the velocity exponent m. The values of other parameters are kept 
fixed at a = 0.7, 6o = 0.002, c — 0.008, and d = 0.0001. As can be verified these equations exhibit a strong volume 
contraction in the four dimensional phase space. We note that there are widely differing time scales corresponding 
to a, bo, c and d (in the decreasing order) in the dynamics of the model. For this reason, the equations are stiff and 
the numerical integration routines were designed specifically to solve this set of equations. We have used a variable 
order Taylor series expansion method as the basic integration technique where the coefficients are determined using 
a recursive algorithm. 



III. SUMMARY OF BIFURCATION EXHIBITED THE MODEL 



The model exhibits a rich variety of dynamics such as period bubbling, period doubling, and complex bifurcation 
sequences referred to as mixed mode oscillations(MMOs) in literature. Here, we will briefly recall only those aspects 
of the bifurcation diagram relevant for the discussion of relaxation oscillations. The gross features of the phase dia- 
gram in the (m, e) plane are shown in Fig. 3. In our discussion, we use e as primary control parameter and m as the 
unfolding parameter. For values of m > ~ 6.8, the equilibrium fixed point of the system of equations, denoted 
by {xo,yo, zo, (po), is stable. Both xq and zq are ~ e/2 and yo and (po are independent of e. At m = m^, we have 
a degenerate Hopf bifurcation as a function of e. For values less than irid, we have a back-to-back Hopf bifurcation, 
the first occurring at e = and the reverse at e = ec2- The periodic orbit connecting these back-to-back Hopf 
bifurcations is referred to as principal periodic orbit. The dynamics of the system is essentially bounded by these 
two Hopf bifurcations. In Fig. 3, the broken line represents the Hopf bifurcation and the dotted lines correspond 
to the locus of first three successive period doubling bifurcations. The inner, continuous lines represent the locus 
of saddle node bifurcations corresponding to period 3, 4 and 5 which are the first three dominant periodic windows 
in the alter nating periodic chaotic bifurcation sequence. Complex bifurcation sequences, characterized by alternate 
periodic-chaotic sequences are seen in the hatched region of the parameter space. A codimension two bifurcation point 
in the form of a cusp at (ec, mc) formed by the merging of the locus of two saddle node bifurcations of the principal 
periodic orbit (represented by bold lines) is shown as filled diamond in Fig. 2. Phase plots and other diagrams have 
been obtained by plotting the maxima of any one of the variables x,y,z or as a function of the control parameters 
(e, to). 



IV. MECHANISM OF RELAXATION OSCILLATIONS 



One characteristic feature of the dynamics of the system is its strong relaxational nature. This feature persists even 
in regions of the (to, e) plane wherein complex periodic-chaotic oscillations are seen (hatched region in Fig. 3). The 
presence of relaxations oscillations and complex periodic chaotic oscillations are interrelated and are a result of the 
geometry of the slow manifold. (For details see Ref. |^^.) Relaxation oscillations that manifest in the model is a type 
of relaxation oscillation wherein the fast variable takes on large values for a short time after which it assumes small 
values of the same order of magnitude as that of the slow variables. The time spent by the fast variable in the part 
of phase space where the amplitude is small is a substantial portion of the period of the orbit. A typical plot of x{t) 
( continuous line) and z{t) ( dotted line) are shown in the inset of Fig. 4 for e = 200.0 and to = 1.2. 
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To understand the nature of the relaxation oscillations, we first study the structure of the slow manifold (5) and 
the behavior of the trajectories visiting different regions of S. The slow manifold of a multiple timescale dynamical 
system is given by the surface spanning the time invariant solutions of the fast variable. In our case, it is given by 

X = g{x, y, (j)) — ~box^ + xS + y — (10) 

with 6 — (f)"^ ~ y — a. Here, the slow variables y and (p (and therefore 5) are regarded as parameters. Further, as we 
will see below, it is simpler to deal with the structure of the slow manifold in terms of the 5 instead of both y and 0. 
Then, the physically allowed solution of the above equation is 



6+^5^+ Abpy 

^0 ^''^ 

where 6 can take on both positive and negative values. Noting that bo is small and therefore ^ 46oy, two distinct 
cases arise corresponding to (5 > and (5 < for which x ^ S/bg and x ^ —y/6 respectively. Further, since the slow 
variable (p and y take on values of the order of unity, the range oi S = S{y, </>) is of the same order as that of 4> and 
y (as is evident from Figs. 4 and 5). Thus, we see that x ~ —y/S is small and x ^ 6 /bo is large. For values around 

1 /2 

S — and positive, we get x ^ (y/bo) 



The bent-slow manifold structure along with the two portions of the slow manifold, namely, Si {S > 0) and 5*2 {S < 0) 
are shown by bold lines in the {x, S) plane in Fig. 4. We have also shown a trajectory corresponding to a mono-periodic 
relaxation oscillation (m = 1.2 and e = 200.0) by a thin line. As can be seen, the trajectory spends most of the time 
on S'l and 5*2. A local stability analysis for points on Si and 5*2 shows that dg/dx = 6 — 2box is negative implying 
that the rate of growth of x is damped. Hence these regions, S'l and 5*2 will be referred to as attracting or "stable". 
For points below the line 2box = 6 {6 > 0), dg/dx > and hence we call this region as 'unstable' (shaded region of 
Fig. 4). Even then, the trajectory starting on 5*2 does continue in the direction of increasing S beyond S — 0. We note 
that this region is not a part of the slow manifold. Once the trajectory is in this region, it moves up rapidly in the x 
direction (due to the 'unstable' nature) until it reaches x = 6/2bo Hue, thereafter, the trajectory quickly settles down 
on to the S'l part of the slow manifold as dg/dx becomes negative. As the trajectory descends on Si approaching 
S2, we see that the trajectory deviates away from S'l. This happens when the value of x is such that 2box < 6, i.e., 
dg/dx > 0. Thus, points on S'l satisfying this condition are locally unstable. Thus, the trajectory makes a jump 
from S'l to S'2 in a short time. This roughly explains the origin of the relaxation oscillation in terms of the reduced 
variables S and x. 



The actual dynamics is in a higher dimensional space and a proper understanding will involve analysis of the movement 
of the trajectory in the appropriate space. Moreover, unlike the standard S— shaped manifold with upper and lower 
attracting pleats with the repulsive (unstable) branch p3|| , in our model, both branches of the bent-slow manifold are 
connected, and there is no repulsive branch of the slow manifold. Thus, the mechanism of jumping of the orbit from 
S2 to Si is not clear. In order to understand this, consider a 3-d plot of the trajectory shown in Fig. 5. The region 
S'2 corresponding to small values of x lies more or less on the y — 4> plane and the region S'l corresponding to large 
values of x is nearly normal to the y — (j) plane due to the large b^^ factor. Regions S'l and S2 are demarcated by the 
'fold curve' given hy 5 = cf/^ — y — a = which dominantly lies in the y — (f> plane. The rapidly growing nature of the 
trajectory lying to right of the 'fold curve' is due to dg/dx > 0. 

The principal features of the relaxation oscillations that we need to explain are: a) very slow time scale for evolution 
on S'2, b) fast transition from S2 to S'l and c) evolution on S'l. In order to understand this, it is necessary to establish 
how the trajectory (viz, x,y,z and (jj) visits various regions of the slow manifold in a sequential way. However, our 
emphasis is more on those aspects of relaxation oscillation pertaining to the issue of time scales raised in the Intro- 
duction, i.e., the timescales involved in the stress - time curve (Fig. 1 ) . (A detailed investigation on the behavior of 
trajectories on this slow manifold has also been carried out. See Ref. [p2[.) However, to understand the dependence 
of stress 0(t), we would also require information of y which in turn depends on z. In order to understand this, we 
shall analyze Eqs. (7) and (9) by recasting them in terms of 6 in various regions of S. This will help to understand 
the general features of the flow viz., on S'2, just outside S2 and on S'l. In the whole analysis, it would be helpful to 
keep in mind the range of values of x,y, z and ^, shown in Figs. 4 and 5, in particular their values as the trajectory 
enters and leaves S'l. 

First, consider rewriting Eq. (7) valid on the slow manifold S in terms of S: 

y — bo{xS — xy + az). (12) 
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The presence of the z variable in Eq. (12) poses some problems. Using detailed arguments based on the knowledge 
of the magnitude of x and z just inside, on and outside 5*2, it can be shown that the trajectory enters 5*2 at small 
values of y and leaves ^2 at relatively larger values. Further one can show that there is a turning point for y on Si 
(see Fig. 5). For details, see Ref ||. 

With this information on the evolution of y on 5*2, we now consider the changes in as the trajectory enters and 
leaves 5*2. From Eq.(9), it is clear that an yield drop starts when x is large (i.e., when x ~ 5/6o on S\) and ends 
when X is close to minimum, when the trajectory is on 5*2, which implies that is small when the trajectory enters 
Si- Using the value of a; = yl\^\ on 5*2 in Eq. (9), we find that e 3> (/'™2//|(5|, since y is near its minimum value as 
the trajectory enters ^2. Thus, increases linearly from small values of at a rate close to de <C 1. We recall that 
the loading rate in the experimental stress-strain curve was Ke'aide in scaled variables), which now can be understood 
as due to the structure of the slow manifold. This is a direct consequence of the fact that the magnitude of x remains 
constant as i: ~ for the entire interval the trajectory on 5*2. This is consistent with what we argued from the stress 
- time plot (Fig. 2), namely, the mobile dislocation density should be constant during the loading period. As the 
trajectory moves into 5*2, y goes through a maximum whereas (j) continues to increase since x ~ y/\S\ remains small. 
However, as the trajectory is just outside 5*2 for which x (y/bo)^^^ for J > and small, 0™ [y/bo)^^'^ ~ e, since 
and y are relatively large which implies that is about to decrease. The above discussion on y and for region just 
outside and inside the fold curve also gives us the direction of movement of the trajectory in this region, namely, it 
enters 5*2 in the region corresponding to small values of y and 0, and makes an exit for relatively larger values of 
and y (compared to their values as the trajectory enters 5*2). Further, as i: 0, we see that the dynamics on ^2 is 
controlled by the slow variables. 

Finally, just to the right oi S = line, x ~ xS, with S very small, which suggests that the time constant is small. 
Thus, the growth of x is slow in the neighbourhood of (5 = 0, and is tangential to the S2 plane even in the 'unstable' 
region. However, once the trajectory moves away from 6 = 0, the growth of the trajectory is controlled by dg/dx 
and hence the time scale of growth of x is of the order of S^^ which is of the order of unity explaining the short time 
span of the stress drop seen in Fig. 2. This also explains why the trajectory tends to leave stable portion of the slow 
manifold ^2 and move into the 'unstable' region. 

Once in the 'unstable' region, the value of x continues to grow in this region of the phase space as can be seen from 
Eq. (9) until the value of x is such that 0™x = e is satisfied. Beyond this value of 0, is negative. Thus, the 
trajectory leaving Si eventually falls onto the Si part of the slow manifold. We can again evaluate y and just as the 
trajectory reaches 6*1. Using x ~ S/bo in Eq. (7), it can be shown that y decreases. Now, consider the equation for 
0. Using a; ~ ^/6o on 5*1, we see that 0'"^/6o > e. Thus < when the trajectory reaches Si with a time constant 
~ d/bo which are relatively fast. (These statements are true only as the trajectory hits 5*1.) We recall here that in 
the experimental time series, the stress drops from a peak value to its minimum in a very short time span. Further, 
we have argued that this should be the sum of contributions arising from fast multiplication of dislocations (which 
we have already argued has a time scale of S~^) and subsequent immobilization. The latter is reflected in another 
rapid time scale ~ d/bo. This explains the difficulty in separating the contributions arising from the two processes in 
the experimental time series. Moreover, since x is a fast variable, the changes in x component dominates the descent 
of the trajectory. Finally, as the trajectory approaches ^2, dg/dx becomes positive and the trajectory jumps from 
Si to 5*2. Combining these results, we see that the trajectory moves towards the region of smaller values of y and 
entering Si in a region of small values of y and 0. 

In summary, the sequential way the orbit visits various parts of the phase space is as follows. The trajectory enters Si 
part of the slow manifold in regions of small y and making an exit along ^2 for relatively large and y. Thereafter, 
the trajectory moves through the 'unstable' part of the phase space before falling onto the 5*1 and quickly descends 
on Si . This completes the cyclic movement of the trajectory and explains the geometrical feature of the trajectory 
shuttling between these two parts of the manifold and the associated time scales. 

Now, the question that remains to be answered is — do the trajectories always visit both and Si or is there a 
possibility that the trajectory remains confined to 5*1 ? It is clear that if the former is true, relaxation oscillations 
with large amplitude will occur and if the latter is true, these are likely to be nearly sinusoidal small amplitude 
oscillations. Here, we recall that the coordinates of the saddle focus fixed point are = zq ~ e/2 which is much 
larger than the value oi x ^ y/ \6\ on Si. Thus, the fixed point located on the Si will be close to the 'fold' at the 
first Hopf bifurcation which occurs at small values of e = '--^ 5. Due to the unstable nature of the fixed point, the 
trajectories spiralling out are forced onto the ^2 part of the manifold resulting in relaxation oscillation. This point 
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has been illustrated by considering the example of a period eleven orbit for m = 1.2 and e = 267.0 shown in Fig. 6. 
As is clear from this diagram, the small amplitude oscillations are located on the Si. As e is further varied, the small 
amplitude oscillations grow with e, but the relaxation nature does not manifest until the orbit crosses over to 82- To 
the best of the authors knowledge the mechanism suggested here for pulsed type relaxation oscillations is new. 

As we will see, the analysis of the slow manifold and the time scales operating in different parts of the phase space 
will be useful in providing an appropriate interpretation of the various branches of the SRS. 
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V. NEGATIVE STRAIN RATE SENSITIVITY 



At the outset, we stress that it has been recognized that the negative unstable branch is not accessible to the dynamics 
of the PLC effect. Even so, early formulations and the way experimental measurements have been carried out has 
given rise to considerable confusion. The purpose of the material presented below is to briefly discuss the concept of 
negative SRS and working methods adopted in the literature, and also clear some misconceptions. 

Theories of dynamic strain ageing (DSA) assume that the interaction of dislocations with solute atoms when averaged 
over the specimen dimensions can be represented by a constitutive relation connecting stress, strain and strain rate 
which is conventionally written as 

cf = he + F{e). (13) 

The basic assumption inherent in Eq. (13) is that stress can be split into a function of e and another of e alone. Then, 
the SRS is defined as 



dlr 



4 (»' 



Clearly, this definition uses e as a state variable. This unfortunately is not correct since strain is history dependent. 
Inspite of this, conventionally, strain is fixed at a small nominal value and the flow stress at that value is used to 
obtain the SRS. It is interesting to note that the existence of critical strain for the onset of the PLC effect implies 
that when the nominal strain value is lower than ec, there are no serrations even when the applied strain rate value is 
in the domain of the PLC effect ( < e < Cca). Yet, the onset of serrations for higher strains is somehow reflected 
in the measured nonmonotonic behavior of the flow stress pl| ]. In experiments, by fixing e at some nominal value less 
than Ec, the flow stress ( at the flxed strain) is found to increase as a function of applied strain rate e for e < e^, 
shows a decreasing trend for < e < and again reverts to an increasing trend for e > Thus, the flow 

stress has the form shown in Fig. 1. No explanation has been offered in the literature as to why this nonmonotonic 
behavior should be seen. However, explanation from the dynamical point of view is fairly straightforward and is as 
follows. We recall here that the model predicts the existence the critical strain Cc, as also the existence of window of 
strain rates Cd < e < ec2 within which serration can occur. Thus, from Eq. (9) it is easy to understand the increasing 
order in which the stress-strain curves are placed for increasing values of e when e < and e > . In this range of 
e, the fixed point is stable and thus all trajectories converge to the fixed point. However, for < e < ec2, we note 
that serrations result only for large enough strains, i.e., once the time of deformation is such that strain crosses In 
our theory, serrations are equated with the existence of periodic (or aperiodic) solutions when < e < ecj. These 
steady state solutions are usually reached only after transients die down. Thus, low value of nominal strain implies 
short evolution time which in turn implies that the stress is being monitored at a transient state. Thus, the decreasing 
trend of the fiow stress for < e < is a reflection of the impending periodic (aperiodic) steady state that will 
be reached eventually. Indeed, this was the method followed in our earlier calculation since the procedure was easy 
to implement numerically However, in many experimental situations, it is not possible to choose a nominal strain 
value low enough that it is less than for the entire range of strain rate values. In such a case, since the stress-strain 
curves are serrated, there is an ambiguity in the value of stress to be used. A working method adopted is to use a 
stress value as the mean value of the upper and lower stress values Then, the flow stress appears to decrease for 
the domain of applied strain rate values where the PLC effect manifests. Thus, this method gives the impression of 
actually measuring the unstable branch. 



The above methods are not suitable for adoption since they do not permit the use of the knowledge of the slow 
manifold. There is an alternate method which uses the relaxation oscillations inherent to the dynamics of the PLC 
effect. In this method, by analogy with electrical analogues, one assumes that there exists a family of curves F{lp) 
for each e of the form shown in Fig. 1 which trigger relaxation oscillations in the form of plastic strain rate bursts and 
stress drops. By comparing the measured stress drops and strain bursts, one concludes the existence of the unstable 
branch, but one never records any points in this region. This method is suitable for our study since we will use 
relaxation oscillations arising in the model. 

In the following we shall argue that the two slow time scales in the dynamics actually translate into the two stable 
dissipative branches of the SRS and the two fast time scales to jumps in plastic strain rate across the stable branches 
of the negative SRS. Since SRS represents (/) as a function of the plastic strain rate kp = (j)"^x, in Fig. 7, we have shown 
a projection of the phase space trajectory on the (f> — ip plane (instead of </i — a;) corresponding to a monoperiodic 
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relaxation oscillation. ( Here, we have retained the same notation for the scaled plastic strain rate as for the unsealed 
one.) The unstable fixed point is also shown. Starting from any initial value around the unstable focus, trajectories 
spirals out converging onto the limit cycle. In Fig. 7, we have identified different regions of the phase space with 
different regions of the slow manifold, and 5*2. We first note that there is a considerable similarity between Fig. 
7 and the schematic representation of the relaxational oscillation obtained using the negative SRS shown in Fig. 1. 
Note also that in contrast to the artificial flat parts B'C and CD' of Fig. 1, the equivalent parts in Fig. 7 have a 
finite negative slope. Lastly, as in experiments, the strain rate jump from B to C is over two orders of magnitude. 

Here, we set up a correspondence between the dynamics in the phase space (Fig. 7) and the slow manifold (Figs. 4 
and 5). From our earlier discussion, we know that when the trajectory is on £'2, x is constant and small in magnitude. 
Consequently, according to Eq. (9), (f) should increases linearly and hence this corresponds the rising branch AB in 
Fig. 7. Further, noting that x Q for the entire interval of time spent by the trajectory on ^2 (sec Fig. 4), the branch 
AB of Fig. 7 corresponds to the pinned state 0/ dislocations. For this branch, one can easily see that the (mean) value 
of <S ~ 3.5 using Eq. (14). Further, as we move up on this branch towards B (Fig. 4), the value of 5 approaches 
zero, and </> reaches its maximum value. Once 5 becomes positive, the trajectory leaves 5*2, and thus, the strain rate 
jump from B to C in Fig. 7, corresponds to the trajectory jumping from S2 to Si in Figs. 4 and 5. Note that the 
slope d(j)/dip for this portion of the orbit is quite small and negative unlike the zero value for the equivalent part in 
Fig. 1. Further, we know from Fig. 4, once the trajectory reaches Si, the value of x decreases rapidly resulting in 
the decrease of e^. Thus, the region CD in Fig. 7 corresponds to the movement of the trajectory on Si ( in Fig. 4 
for which 5 > 0). For this branch, one can quickly check that the strain rate sensitivity S is positive having a mean 
value ( ~ 1.5) which is a factor of two less than that for the branch AB, implying that the nature of dissipation is 
quite different from that operating on AB. This is consistent with known facts about the two branches as mentioned 
in the introduction. Combining this with the fact that x is decreasing, the branch CD in Fig. 7 mimics the equivalent 
branch CD' in Fig. 1, which is identified with the slowing down of the mobile dislocations without solute atmosphere. 

We recall that the stress drop duration has contributions from two fast processes, namely, dislocation multiplication 
and it's subsequent immobilization. But, these two time scales could not be separated in the stress strain curve. How- 
ever, in the present phase plot representation (Fig. 7), we sec that the fast multiplication of dislocations correspond 
to BC and that of immobilization to CD. This correspondence has been possible due to the mapping of the relevant 
time scales in the dynamics of the dislocations obtained from the analysis of the slow manifold to the various regions 
in the phase plot thereby allowing us to identify the individual contributions. (Note also that in Fig. 7, we have 
plotted points of the trajectory at equal intervals of time which shows that the time interval corresponding to BD is 
small.) Prom Figs. 4 and 5, we see that, as the trajectory descends on Si part of the slow manifold and gets close to 
^2, it leaves Si, since dg/dx becomes positive {x 50). Further, the strain rate sensitivity parameter S changes sign 
at D. For the corresponding DA part in Fig. 7, the slope is small and negative as for the part BC. Noting that B and 
D are the points at which S turns negative, and noting that the fixed point is unstable, the so called "unstable branch" 
of the SRS, not accessible to the dynamics, can be inferred by drawing a (dotted) line connecting the maximum and 
the minimum of the stress and passing it through the unstable fixed point (Fig. 7). 

We will now attempt to use the results of our analysis of time scales in different regions of the slow manifold to obtain 
the dependence of Cp as a function of (j). The equation for Cp is 

^ =xm(f)"'-'^ 4) + X (15) 

which on using Eq.(6) and (9) gives 

dip ^ ep{^+S)-elV-f + ^) + y,^"^ 

d<P d{e-ip) ^ ' 

Here we note that in the slow manifold description, all slow variables appear as parameters. However, since SRS 
describes the dependence of the slow variable (j) as & function of the (derived) fast variable Cp, we will consider the 
other two variables y or 5 or both as parameters. Numerical solution of Eq.(16) has been attempted using y and 
5 as parameters. Good numerical approximation is obtained by noting that y and hence 5 is periodic. Thus, any 
reasonable approximation for the periodicity of y, for example, sine function with a proper amplitude and phase gives 
a good fit with the phase plot. However, our interest here is to obtain approximate expressions for ep{(j)) on different 
branches. For this reason, we will use typical values of 5 and y for the interval under question. Prom section 4, the 
trajectory has different dynamics in different regions of the slow manifold. These are : (1) on ^2 where x is nearly 
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zero for the entire time spent by the trajectory on 5*2, (2) just outside S2 where x ~ a;(5,(3) on Si where x ~ S/b^ for 
e'p > e and (4) when the trajectory jumps from Si to S'2. Approximate solutions obtained for these cases are shown in 
the phase plot by solid lines. Details are given in the Appendix. It is clear that these solid lines reproduce the general 
features of the phase plot quite well. We stress here that these lines correspond to the simplest approximation. 
The above analysis refers to a fixed value of e. As a function of e, we find that the magnitude of the stress drops 
increases initially and then decreases. This feature is a direct result of the existence of back-to-back Hopf bifurcations 
in the model. On the other hand, experimentally one sees only a decreasing trend. While the decreasing trend is 
consistent, the increasing trend seen in the model for low strain values can be traced to the effect of another crucial 
parameter in the model, namely b^. We recall that this parameter corresponds to the remobilization of immobile 
dislocations. For the value of &o used in the present calculation, the bifurcation from the steady state is a mildly 
subcritical Hopf bifurcation, i.e., across the transition the amplitude of the stress change is abrupt but the magnitude 
is small. However, for smaller values of 60 1 this jump can be made sufficiently large in which case the amplitude of 
the stress drops can be made to decrease with e right from the onset of the PLC effect. 



VI. DISCUSSION AND CONCLUSIONS 

The study of the relaxation oscillations in the model was motivated by the need to explain the apparent inconsis- 
tency between the time scales observed in experimental stress-time series and those that could be argued on the basis 
of the negative SRS feature commonly used in the literature. The study of relaxation oscillations using the geometry 
of the slow manifold has helped us to identify different time scales operating in different regions of the phase space, 
apart from showing that the nature of the relaxation in the model is due to the atypical bent geometry of the slow 
manifold. This geometry is very different from the standard S— shaped manifold and hence the relaxation oscillations 
seen here differs qualitatively from that seen in systems with 5'— shaped slow manifold. Some comparative comments 
between these two types of manifolds may be in order here. As in the 5— shaped manifold, there are two attracting 
branches in our case also, namely Si and S'2. The dynamics on S'2 is slow as it is controlled by the slow variables y 
and (f>. On the other hand, on Si, the time dependence of the trajectory is largely controlled by the fast variable x. 
In this sense, the dynamics on S2 is slow and that on Si is fast, with a large magnitude and with a much smaller 
magnitude from Si to S2. Though there are two fast jumps as in the S— shaped manifold, in our case, there is no 
equivalent unstable part of the slow manifold which causes these jumps. 

The analysis of the time scales controlling the relaxation oscillations has been directly used to reconstruct the relax- 
ation oscillation in the (t> — ip plane which bears a strong resemblance to the the relaxation oscillations resulting from 
the assumed form of the negative SRS. The information on different time scales operating in different regions of the 
slow manifold has been used to calculate the dependence of (j) on ip for the two dissipative branches and the associated 
strain rate jumps between them. This has helped to identify the various regions of the slow manifold with the stick 
state and the slip state of dislocations. It has also helped us to clarify the inconsistency in the time scales of the 
dynamics. Further, several important features of the SRS derived from the model compare well with that reported 
in the literature. In particular, we note that the slope of the first dissipative branch AB is larger than that of the 
second branch CD (Fig. 7). Further, we recall that —5 = —cjf^ +y-\-a which is positive for AB, gradually approaches 
zero as B is reached followed by strain rate jump. Similarly, for the branch CD, 5 approaches zero as we approach D 
followed by a jump in the strain rate. Thus, vanishing of 5 is indicative of strain rate jumps just as the strain rate 
sensitive parameter S also vanishes. Noting that y is the immobile dislocation density, it is tempting to interpret 5 

1 /2 

as being related to some kind of effective stress. (Recall that the effective stress is a* ~ cFa 

- HN-^, where H IS 

the work hardening coefficient.) Thus, the points at which strain rate jumps occur correspond to points at which the 
effective stress vanishes which is very much like the classical explanation. Since the definition of strain rate sensitivity 
assumes strain as a state variable which is not true, 5 may be an effective alternate parameter for defining strain rate 
sensitivity. Thus, it is nice to see that we can attribute a physical meaning to this parameter. 

The analysis has also helped us to provide a dynamical interpretation of the negative SRS. The analysis also shows 
that the large jumps in the strain rate across the stable branches are due to the relaxational nature of the dynamics 
which in turn is a result of bent nature of the slow manifold and the fact that the bifurcation is of the Hopf type. 
Using this, we have inferred the existence of the unstable branch as containing the two points, B and D, where strain 
rate jumps (where 5 and S are zero) and the unstable fixed point. In this sense, Hopf bifurcation is at the root of the 
'negative' SRS. Similar features of SRS were found to operate in a model designed to mimic stick-slip dynamics of 
tectonic faults |33|. There are other studies on stick-slip dynamics, both experimental [p6| and theoretical p^, which 
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support the view that Hopf bifurcation was found to be responsible for the instabihty. Thus, it is hkely that Hopf 
bifurcation is relevant to situations where stick-slip dynamics operates and wherever one measures the two stable 
branches and the jumps across the branches | |3^ . 

The relaxation oscillations in the model are reminiscent of the canard type of oscillations in multiple type scale dy- 
namical systems [ ^9|j4C[| . The latter type of oscillations result from 'sticking' of the trajectory to the repelling part of 
the 5— shaped slow manifold before jumping to the attracting pleat of the slow manifold. In our case, although these 
oscillations have a similarity with canard type of solutions, the repelling part of slow manifold does not exist. Instead, 
the trajectories stick to the 'unstable' part of the phase space where the dynamics is accelerated once the trajectory 
moves well into this region. This aspect coupled to the fact that there is no inherent constraint on the manifold 
structure leads to oscillations of all sizes. It is clear that such oscillations result from the trajectory 'sticking' to the 
direction of the 5*2 plane and moving into the 'unstable' part of the phase space by varying amounts each time the 
trajectory visits 82- These jumps translate into stress drops of varying sizes which are generally seen in experimental 
time series (Fig. 2). This also means that e'p — (f> is not a simple limit cycle, and the simplistic approach of inferring 
the 'negative' SRS should be given up. The present analysis stresses the importance of using sound dynamical tools 
such as the slow manifold as the basis for studying more complex oscillations rather than phenomenological concepts 
such as the negative SRS. 



APPENDIX A: 



Here, we obtain approximate analytical expressions for ep{(j)) for different regions of (/> — e'p phase plot (Fig. 7) using 
the knowledge of time scales obtained from the analysis of relaxation oscillations. For the numerical evaluation, the 
values of control parameters have been chosen as e = 200, m = 1.2, bg = 0.002 and d = 0.0001. 

Region AB : When the trajectory is on S2, i ~ 0, for the entire interval of time. Using x = — 2//<5 in Eq.(15), we get 

dip _ (j)'''-^my 

Noting that S = 0™ — y — a, this equation can be integrated thereby reducing the number of parameters to one, 
namely y. Integrating, we get 

where ep(0) and (j){0) refer to their respective values as the trajectory enters 82- In Fig. 7, we have used ep(0) = 4.7, 
0(0) = 2.2 and y = 6.15. 

Region BC : This region corresponds to the jump from S2 to Si. This happens when the trajectory is just outside 

1/2 

82- For this region, is near (pmax, and since x ~ (y/^o) j the evolution of x is well described by i '--^ x5, implying 
that the time of evolution is very short. Thus, we can regard the evolution of <f> as being mainly determined by that 
of x. (This region also corresponds to (5 > and small ^ 0.2.) Thus, we use — (pmax on the RHS of Eq. (9). Then, 

^ = d(e-Ca.^..e**) (A3) 

where Xs2 is the value of x at the time of leaving 82- Integrating this equation with initial conditions at t — 0, 
(j) = 0(0) = (t>max, we get 

git ^ _ {(f>{t) - (f>rnax)5 ^ 

T fh™ T rh"^ d 

•^S2Vmax ■'^s-l^max^ 

Clearly, the first term is small since the time span of evolution that we are interested is also ~ S. Now consider Eq. 
(15). Using X ~ xS, we get 



dcp , ( med \ mf^d 
—£- = e ' ' ' 

dt 



(A5) 
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Since med/<j)max « S, we drop the first term. Integrating the above equation with the initial value ep(0) = e, leads 
to 

max 

^ m<P(l)'^^^Xs^{(t)max5/tnd - e) - - (j)max)emd 
In Fig. 7, we have used the values 5 = 0.021, (t)max = 4.98 and = 1-7. 

Region CD : Consider the trajectory on with x ~ 6 /bo and ip ^ e. Then, Eq. (16) reads, 

Since e/lp < 1.0, wc expand the denominator and retain terms upto e^^ . We note hero that on 5*1, x is rapidly 
decreasing and therefore, using slow manifold values is not a good approximation. Even so, as a simplest approximation 
we use X ~ 5/hQ. Then, we get 

dip ^ me^ho eb^y m4r--'^5 bpy 

d(f) ~ (^""+^6 (p'^dS^ bo d6 ^ ' 

In this equation, both y and 5 appear as parameters whose values are chosen appropriate to this region. Integrating 
the above equation with the initial values of 0(0) and ep(0), we get 



(m - l)rf^2 



+_(^™ - 0^(0)) - -^{cp - m)) + ep(0) (A9) 



In Fig. 7, we have used the values of S = 3.57, y = 1.0 with (/)(0) = 4.3 and ep(0) = 7460.0. 

Region DA : This region again corresponds to the the trajectory descending on but ip < e. Then, Eq.(16) can be 
written in the form 

de, ^ ep{^+5)-ep'{^ + ^)+y<P"^ 
# e(i(l-^) 

In this region, the value of cp is slowly varying with its value near the minimum for which 6 ~ —0.15. Since <p is near 
^min, we use (p = 4>min and regard the variation as largely arising due to the changes in Cp. Using (1 — ^)~^ (1 + ^), 
we have 

^ = -{A^ip + B.cp^ - Cie-/) (All) 

where = ^ - ^ and Bi = - and Ci = -f^ + ^2^$^- Since Ci < Ai and Bi, we drop the last 

term. Integrating the above equation with the initial conditions, 0(0) = (pmin and e'p = e, we get 

= Ai+eBi(l-e^i(^-^™.")) ^"^^^^ 

In Fig. 7, we have used 4>min = 1-55. 
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FIG. 1. Schematic plot of the SRS. Branch B' D' shown by the dashed Une describes the negative strain rate sensitivity of 
the PLC effect. Dotted hnes represent the discontinuous strain rate jumps leading to serrations in the stress strain curve. 
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FIG. 3. Phase diagram of the model in (m, e) plane. See text for details. The broken line corresponds to the locus of Hopf 
bifurcations and dotted lines to the successive period doubling bifurcations. See text for details. 



16 




FIG. 4. Evolution of a trajectory (thin lines) along with the bent-slow manifold {Si and S2 shown by thick lines) structure 
in the x ~ 5 plane, for m = 1.2 and e = 200. Inset shows the time series of the x variable (continuous line) and z variable 
(dotted line). 
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FIG. 5. Evolution of the trajectory along with bent-slow manifold {Si and S2) structure in {x,y, cj)) space indicated by the 
gray plane, for m = 1.2 and e = 200.0. 
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FIG. 6. Evolution of the trajectory along the bent-slow manifold (Si and 52) structure for m = 1.2 and e = 267.0. 
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FIG. 7. Empty circles show the phase space projection of a vs ep corresponding to a relaxation oscillation. The unstable 
fixed point is shown by a filled diamond. The dotted line through the fixed point represent the apparent negative SRS region. 
The thick lines are analytical approximations of corresponding regions. 
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